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ABSTRACT 


A localized perturbation of a magnetic flux tube produces a pair of wave trains 
that propagate in opposite directions along the tube. These wave packets disperse 
as they propagate, where the extent of dispersion depends on the physical properties 
of the magnetic structure, on the len gth of the in i tial e xcitation, and on its nature 
(e.g., transverse or axisymmetric). In lOliver et al.l (120141 1 we considered a transverse 
initial perturbation, whereas the temporal evolution of an axisymmetric one is examined 
here. In both papers we use a method based on Fourier integrals to solve the initial 
value problem. Previous studies on wave propagation in magnetic wave guides have 
emphasized that the wave train dispersion is influenced by the particular dependence of 
the group velocity on the longitudinal wavenumber. Here we also find that long initial 
perturbations result in low amplitude wave packets and that large values of the magnetic 
tube to environment density ratio yield longer wave trains. To test the detectability 
of propagating transverse or axisymmetric wave packets in magnetic tubes of the solar 
atmosphere (e.g., coronal loops, spicules, or prominence threads) a forward modelling 
of the perturbations must be carried out. This is left for a future work. 


Subject headings: Sun: corona - Sun: magnetic fields - Sun: oscillations 
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INTRODUCTION 


The solar corona supports a rich variety of wave phenomena whose main features are deter¬ 
mined by the magnetic field and plasma structuring. For our purposes it is important to distinguish 
between different wave manifestations such as driven and impulsively excited perturbations. Driven 
waves are often caused by a continuous, periodic exciter at the photospheric or chromospheric level 
and propagate away from the place where they are generated. This is the case, for example, of 


waves propagating upwards along sp icules (lOkamoto &: De Pontieu 


or along the coronal magnetic field (jTomczvk et al 


2007 


2011 


De Pontieu et al.ll2012ll 


Tomczvk fc McIntosh 2000 ). The dis 


turbances propagating in filament threads observed by Lin et al. ( 2007 i Hokamoto et al (2007) 


are 

probably representative of driven waves too, although th e exciting agent is not clear in this case . 
Rapidly propaga-ting quasi-periodic coronal disturbanc es (jLiu et al.l l201ll . l2012l : IShen &: Liul l2012l : 
Yuan et al.l l2013l : IShen et al.l 120131 : iNistico et al.l l2014l l also seem to be caused by a continuous 
driver that in this case i s assoc i ated to a flare or the initiation of a CME. For a modelling of these 
events see lOfman et al.l (|201lh : iNistico et al.l (j201J). When a magnetic flux tube (e.g., spicule, 
coronal loop, or prominence thread) is subject to a periodic driving it shows a distinct response. 
An axisymmetric driver that periodically compresses and expands the tube at one position leads 
to the propagation of a sausage-type perturbation (see Figure [11(a)). A periodic transverse shaking 
of the tube results in a ki nk-type perturbat ion by which the tube boundary acquires the shape of a 
snake (see Figure 1(a) of [Oliver et al.l 120141 . hereafter Paper I). A torsional driver leads to Alfv& 
wave propagation along the tube with no apparent modification of its boundary. 



(a) (b) 

Fig. 1.— Shape of the boundary of a uniform cylindrical magnetic flux tube subject to (a) a fast 
sausage normal mode perturbation and (b) a localized axisymmetric perturbation. 


An impulsive energy release can lead to both standing and propa g ating disturbance s. The 
transverse coronal loop oscillations studied by lAschwanden et al.l (Il999l ) ; iNakariakov et al.l (jl999l ) 
constitute a good example of the first type of phenomenon, for which there are many observational 
studies in the literature. For recent observations using the Atmospheric Imaging Assembly (AIA) 
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on board the Solar Dynamics Observatory, see iNistico et alJ (j2013l h lAnfinogentov et alJ (|2013l ) 


Our interest in this work is in the second type of event, i.e., a magnetic flux tube perturbed by 
a sudden, localized expansion or compression such as that of Figure [T]^b) (see also Figure 1(b) of 
Paper I for a localized transverse excitation). In a dispersive medium, such initial pulse disperses 
and transforms into a wave train that induces quasi-periodic plasma oscillations with a range of 
periods (rather than the single-period oscillations of driven propagating waves or standing oscilla¬ 
tions). Hence, the time signal (e.g., Doppler shift or spectral line intensi ty) collected at a given 
position shows a charact eristic “tadpole” shape in its wavelet diagram le.g.. lNakariakov et al.ll2004l : 


Jelmek &: Karhckvll201Cl b The important difference between driven and impulsively excited waves. 


then, is that the first ones have the period of the driving agent, whereas the second ones display 
a continuous variation of the period that depends on the size of the initial perturbation and the 
physical properties of the medium. 


Several observations of q uasi-periodic oscillations hav e been ascribed to the disp ersive behavior 
of propagating wave packets. IWilliams et all ( 20011 . 120021 1 uKatsivannis et al.l ( 20031 ) used the Solar 
Eclipse Corona Imaging System (SECIS) to record white light and coronal green line images of the 
solar corona during a solar eclipse. They studied the intensity variations in the vicinity of an active 
region and detected periodic variations with period in the range 4-7 s that propagate through the 
loop apex with a speed of 2100 km s“^. The wavelet spectra of signals at different points along 
the wave path display the “tadp ole” shape mentioned above. This wavelet feature i s also present 


in solar decimetric radio bursts (jMeszarosova et al 


2009, 


2011 


Karlickv et al.ll2013h and in radio 


sources lying in a magnetic fan structure (jMeszarosova et al.ll2013l ). 


The impulsive excitation of fast magnetohydrodynamic waves, and their subsequent evo- 


sausage wave trains in a 

coronal slab (h 

lurawski & Roberts 1993a 

d: Nakariakov et al. 2004 

2005: 

Jelmek & Karlickv 

2010: 

Nakariakov et al.ll2012 

Meszarosova et al. 

20141. in a current sheet ( 

Jelmek & Karlickv 

2012:ljelmek et al. 

2012 

Karlickv et al 

2013: 

VIeszarosova et al. 2 

014), and in a magnetic 

funnel 

holes 

(Pascoe et al. 2013 

). Impulsively generated fast sausage and kink wave packets in coronal 


have also been studied bv iPascoe et al.l (|2014l ). Our paper is concerned with the temporal evolution 
of a concentrated axisymmetric compression or expansion of a uniform, cylindrical magnetic tube; 
see Figure [11(b). This initial disturbance leads to the propagation of a fast sausage wave train along 
the tube. The propagation and dispersion of wave packets in magnetic slabs has already been 
investigated numerically in slab geometry, but not in cylindrical geometry. 


This is the outline of our paper. Section [2] contains the governing equations of linear fast 
perturbations in a plasma cylinder in the (3 = 0 limit. In Sections [3] and [4] we give a description of 
fast sausage proper and improper modes, whose spectra are discrete and continuous, respectively. 
The initial perturbation is used in Section [5] to obtain analytical expressions for the time evolution 
of the perturbed variables and the results are presented in Section |6l 


































































































- 4 - 


2. EQUILIBRIUM AND GOVERNING EQUATIONS FOR LINEAR 
PERTURBATIONS IN A COLD PLASMA CYLINDER 

We consider a homogeneous cylindrical magnetic tube embedded in a nniform environment. In 
this work we use cylindrical coordinates (r, (/?, z), with the z-axis pointing along the cylinder axis. 
The magnetic field is uniform in the whole space and parallel to the z-axis, and so can be written 
as B = Bez- The equilibrium density is 


Po{r) = 


Pi, r <a, 
Pe, r > a, 


( 1 ) 


where the subscripts ‘f’ and ‘e’ denote internal and external values, respectively, and a is the tribe 
radius. 

Next we consider small perturbations about this equilibrium, governed by the lin ear, ideal mag- _ 

netoh ydrodynamic equations, in the cold plasma approximation (Equations (l)-(3) of iRuderman &: Erdelvi 
20091 1. In our work v is the velocity perturbation, ^ the Lagrangian plasma displacement, b the 
magnetic field perturbation, P = Bbz/po the magnetic pressure perturbation (with /tq the mag- 
netic permeability of free s pace), and p the density perturbation. Using Equations (1)"(3) of 
Ruderman fc Erdelvil ([20091) it is straightforward to see that, if ^z and Vz are initially zero, then 
they vanish for all times. In addition, perturbations are governed by 


dP 




--VuP + 4f|> 

Po dz^ 


oz 


( 2 ) 

( 3 ) 

( 4 ) 


P = -pov\V • 


( 5 ) 


V = 


di 

dt' 


( 6 ) 


where we have written = b — bzEz and 


^ d 

V_L = V - e^ —. 
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Moreover, va is the Alfven speed, whose values inside and outside the magnetic cylinder are re¬ 
spectively given by 


4i = 


52 

^Ae = 




POPe 


( 7 ) 


In this paper we consider a magnetic tube that is denser than the environment, i.e., pi > pe and 
VAi < VAe- 


3. PROPER (OR DISCRETE) FAST SAUSAGE MODES 
3.1. Eigenfunctions 

For sausage perturbations, P, p, b, and v have no azimuthal dependence and as a conse¬ 
quence the (/^-components of Equations (l3|), ([1]), and (l6|) dictate that b^, and vanish (provided 
that they are initially zero). In addition, the variables are of the form 


P{t, r, z) = P{r) exp[i(—-|- kz)]^ (8) 

and similarly for p, br, and Vr- The eigenfunctions of fast magnetohydrodynamic (MHD) modes 
in a magnetic cylinder have the azimuthal dependence exp{im(p), with m = 0 for sausage modes. 

Equation ([2]) gives the radial dependence of the density perturbation 


Equation ([3]) leads to 



(9) 


C - 1 

PQv\k^ dr ’ 

where the radial wavenumber, kr, is given by 


( 10 ) 


k 


2 

r 


— k‘^v\ 


( 11 ) 


Equations (jlj) and ([6]) give the radial components of the magnetic field and velocity perturbations 
in terms of the radial plasma displacement. 


br = ikBp, 


( 12 ) 
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Vr = —iuj^r- 

Finally, from Equations ([5]) and (fT0|) we obtain 


(13) 


dr‘^ 


+ -^ + k‘^P = 0. 

r dr 


(14) 


To determine the solution to Equation (jl4h we need to take into account the possible signs 
of inside and outside the magnetic cylinder (the internal and external radial wavenumbers are 
here denoted as ki and ke, respectively). For proper fast sausage modes kvAi < < kv^e and this 

implies that A:? > 0, fcg < 0. Then, the solution to Equation (fHl) regular at r = 0 and that does 
not diverge as r —>■ oo is 


P{r) 


CiJo{kir), r < a, 
CeKo{Ker), r> a, 


(15) 


where Ci and Cg are two constants, Jm is the Bessel function of the first kind and order m, and 
Km is the modified Bessel function of the second kind and order m. Moreover, 


kf = 


LO 


— k'^'i 


Ai 


^Ai 


(16) 


and Kg = — fcg is 


, ,2 _ u2 2 
.^2 _ ^ 

^Ae 

It is important to mention that when using these equations to determine the radial wavenumbers, 
the positive value is taken, i.e., ki > 0, Kg > 0. An algebraic expression for Ci and Cg is obtained 
by imposing the continuity of P{r) at r = a. The remaining constant can be imposed arbitrarily. 
Here we use expressions analogous to those of Paper I, that is 



Ci — (Kgfl), Cg — PeV A^JoikiCl) ■ (1^) 

Now that the function P{r) is determined, all other perturbations can be obtained from Equa¬ 
tions m, cni), m, and (I13p . We write the expression of the radial plasma displacement because 
it will be used later, 
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&ir) = 


k- Ko{Kea)jQ{kir), r < a, 
-K~^Jo{kia)KQ{Ker), r > a, 


(19) 


wher e a prime indicates a derivative. From now on we use the relations (jAbramowitz fc Stegun 
1964l i 


Joix) =-Ji{x), Kq{x) =-Ki{x), 


( 20 ) 


and write Equation (fT9|) as 


Cr{r) = 


-A:- Ko{Kea)Ji{kir), r < a, 
Jo{kia)Ki{Ker), r > a. 


( 21 ) 


3.2. Dispersion Relation, Phase Speed, and Gronp Velocity 

Imposing the continuity of the radial displacement at r = a leads to the dispersion relation for 
proper fast sausage modes, 


Mha) KijKea) ^ p ^22) 

kiJo{kia) KeKo^Kea) 

In this work the dispersion relation is used to obtain the frequency of a proper mode with a given 
longitudinal wavenumber, k. Therefore, we often write the frequency as co{k). 

Moreover, the phase speed and group velocity are defined in the usual manner. 


Cph(^) 


OJ 

V 



(23) 


Although the phase speed is not used in our calculations nor in the interpretation of results, it is 
also given here because it is an important parameter of proper modes. The main features of cu, Cph, 
and Cn of proper sausage modes in a rnagnet i c cylinder h ave been discussed in previous works (e.g.. 


Edwin fc: Robertslll983l : [Roberts et al.lll984l : ICallvlll986l l. They are summarised below. 


Eirst, we discuss some properties of the eigensolutions for the possible sign combinations of k 
and OJ] they will be useful when the method of Eourier integrals is applied to the wave propagation 
problem. It is clear that the two pairs of values {k,oj) and {—k,—u}) correspond to the same 
eigenmode, propagating in a given direction along the 2 ;-axis. It is obvious that the solutions 
presented above have this property: if the signs of both k and u are simultaneously changed. 

















then Equations (USD and (HD show that the internal and external radial wavenumbers are not 
modified. As a result, these new values of the longitudinal wavenumber and frequency also satisfy 
the dispersion relation. In other words, if uj{k) is the frequency of a proper mode, then —oj{—k) 
is also a frequency of a proper mode. Using Equations dSD, ([HD, (USD, (USD, and dUD , we see 
that the perturbations of both modes are identical. Einally, we also have Cph(A:) = Cph(—fc) and 
Cg{k) = Cg{—k). Next let us consider a fixed k: if u}{k) is a solution to the dispersion relation, then 
—io{k) corresponds to the eigenmode propagating in the opposite direction along the magnetic 
tube. Again, the two eigenmodes have the same eigenfunctions, but opposite phase speed and 
group velocity. Thus, in the rest of this section it is enough to concentrate in /c > 0, a; > 0. Both 
positive and negative values of k and u will be taken into account when deriving the solution to 
the initial value problem by the method of Fourier integrals (Section [5D. 


The magnetic cylinder supports an infinite number of fast sausage radial harmonics, charac¬ 
terised by the number of nodes of P{r). All radial harmonics have a node at the tube axis, r = 0. 
This is the only node for the fundamental sausage eigenmode. The first overtone has two nodes 
in the radial direction, etc., and all these nodes are inside the magnetic cylinder. An additional 
feature of proper fast sausage modes i s that they only exist for longitu dinal wavenumbers exceeding 
the cut-off value (jRoberts et al.lll984l : iNakariakov fc Verwichtel l2005l i . 




jojVAi 



(24) 


where joj is the jth zero (in increasing order) of Jo{x). Here j = 1,2,... correspond to the 
fundamental fast sausage mode, its first radial overtone, etc. 

In Figure [2] the frequency, phase speed, and group velocity of the fundamental fast sausage 
eigenmode and its hrst overtones are shown for two values of the density ratio, PijPe- The dispersion 
diagram (Figure[2Da)) is rarely found in the literature. It shows that all proper modes, including the 
fundamental one, have a wavenumber cut-off and that their frequency satishes kvAi < oo < kv\e- 
One can see that increasing the density ratio raises the upper frequency limit kv\e and so more 
overtones are present in this diagram. It is worth to mention that linear waves with kvAe < ^ also 
exist, but instead of discrete values of u these solutions can have any frequency (see Section ED. 
Figure [2Db) has been presented many times in other works and so is not discussed here. It is 
enough to note that the phase speed lies between the internal and external Alfven speeds. Finally, 
Figure EDc) shows the group velocity as a function of k. When the wavenumber approaches its 
cut-off value from the right, Cg rapidly increases and tends to VAe, i-e-j it tends to 3 and 2 for 
the red and black lines, respectively. For larger wavenumbers, the group velocity decreases to its 
minimum and finally increases slowly towards vai- In Appendix!^ we show that the phase speed 
tends to vai from above as ka —)• oo. Moreover, the group velocity tends to VAi from below for very 
large wavenumbers. This behavior of Cph and Cg is in agreement with Figures [2](b) and (c). 
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4. IMPROPER (OR CONTINUOUS) FAST SAUSAGE MODES 

When a fixed k is considered, proper modes have discrete frequencies in the range {—kvAe, kvAe) 
(see purple dashed lines of Figure [2]^ a)). Improper mode eigenfunctions can also be expressed in the 
form of Equation ([8]), but they can have any frequency in the ranges (—oo, —kvAe] and [kvAe^oo). 
Paper I contains a detailed discussion of kink improper modes properties. All the arguments 
presented in Paper I are valid for improper sausage modes and so are not repeated here. The total 
pressure and radial displacement perturbations are given by 


P{r) = 


PivX-ykiJo{kir), 


r < a, 


Pev\^a?kl[CjJo{ker) + CYYQ{ker)], r > a, 


irir) = 


- a^k^k- ^Ji{kir), 


r < a, 


-o?ke[CjJi{ker) + CyYi^ker)], r > a, 


(25) 


(26) 


where Y^ is the Bessel function of the second kind and order m, and Cj and Cy are two constants, 


Cj = 

Cy = 


nake 

2ki 

nake 

2ki 


[-kiJoikia)Yi{kea) + keJiikia)Yo{kea)] , 
[ keJ\{kia')jQ{keY) kiJQikia') J\{keay\ . 


(27) 


All other perturbations can be derived from the expressions of Section 13.11 


5. SOLUTION BY FOURIER INTEGRALS 
5.1. Solution to Initial Value Problem 

In this section we derive the expressions governing the propagation of an axisymmetric impul¬ 
sive excitation such as that of Figure [T][^a) . The procedure to obtain the evolution of the pertur¬ 
bations after the magnetic tube suffers a localized transverse displacement is described in detail in 
Paper I, where all the relevant information can be found; here only the main equations are given. 
The initial value problem is specified by two initial conditions on the radial displacement. 


^r{t = 0,r,z) = f{r,z), 

= 0,r,z) = g{r,z). 


(28) 


To impose these initial conditions we write the Fourier transform of ^r{t, r, z) with respect to z. 
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'OO 


ir{t,r,k) = / Crit,r,z) 


dz. 


(29) 


— OO 


For a given longitudinal wavenumber, ^rit, r, k) can be expressed as a linear combination of proper 
and improper eigenfunctions, 


N 


Ut,r,k) = Y, i,ir,k) 



(30) 


where ojj(k) > 0 and cu > 0 are assumed with no loss of generality. In this expression, the funda¬ 
mental fast sausage mode (j = 1) and its radial overtones (j = 2,3,...) have frequencies ±a;j(A:), 
which are the solutions to the dispersion relation (Equation ()22p l. and their radial displacement, 
^j(r, k), is given by Equation (|2T]) . The contribution of each of these proper modes to Equation (l30]l 
depends on the amplitude A^{k), where the -|- and — signs correspond to propagation in opposite 
directions along the magnetic tube. Moreover, N is the number of proper sausage modes for the 
selected longitudinal wavenumber. For |A;| < kd no proper eigenmodes exist and so = 0, which 
means that the sum is not present in Equation (j30l) : for kd < \k\ < kc 2 only the fundamental fast 
sausage mode exists and then = 1; and so on. Regarding improper modes, their contribution to 
Equation (1301) has the form of an integral because they have a continuous variation of frequency 
for a given k. The improper modes eigenfunction ^i_j(r,k) has been presented in Equation (|26p . 

We next define the Hilbert space of functions on the interval [0, oo) by the scalar product 



(31) 


With this definition, the eigenfunctions satisfy the orthogonality conditions 


( 4 -, 6)=0 

{ij,L) = 0, 

{L,L') = q{u})6{u) - u}') 


(32) 


where io,io' > 0 and 



00 


( 33 ) 
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see Appendix iBl 

To determine the functions A^{k) (j = 1,...,A^) and A^{k) we use the initial conditions. 
Taking f = 0 in Equation (1301) we obtain 


N 

fir, k) = Y^ [a+ (k) + Aj (k) 4(r, k) 




POO 

+ [A^ik) + A-{k)]i^{r,k)du;, 

J\k\vAe 


(34) 


f{r,k) being the Fourier transform with respect to z of f{r,z). Now, differentiating Equation (I30p 
with respect to t and substituting t = 0 yields 


gir,k) 


N 


i=i 


-i^uijik) AUk)-AAk) ^jir,k) 



Auj{k)] Lir,k)duj, 


(35) 


with 5 (r, k) the Fourier transform of gir, z) with respect to 

The next step is to take the scalar product of Equations (|M1) and ([35]) with (,j{r,k). Using 
Equation (1321) . this results in two algebraic expressions for A^ik) that lead to 


with 



A^jk) 

2ojjik)Vj{k) ’ 


(j = l,2,...,iV), 


(36) 


Mtik) = - 


PiKoiKeO) 


'J 

PeJoikiOi) 

I o 


iujjf ±iAkirJiikir) 


10 
roo r 


dr 


POO 

y <' 

Ja 


ujjf ±iAK-erKiiKer) dr, 


(37) 


V,ik) = 


+ 


PiKliKeO) 

kf 

PeJ^iha) 


{kirfjlikir) 


dr 


>0 

{K^rfKU^er)-. 


(38) 
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Analytical ex pressions for the integrals in this formula can be obtained. Using Equation (1201) and 
the relations (jAbramowitz &: Stegunlll964l l 


we obtain 


A. 

dx 

d 

dx 


xJ{{x) = xJq{x) — Ji(x), 
xK[{x) = —xKq{x) — Ki{x), 

[x^ [Joix] + Ji (x)] } = 2xJ^{x), 
[K^{x) - Kf{x)] } = 2xKl{x). 


With the aid of these formulae we get 


(39) 


(40) 


f 


rir 

{kirf Jl{kir)— = [Jl{kia) + Jf (fejo)] 
tZ 

- kiaJo{kia)Ji{kia), 


(41) 


f {Kerf Kl{Ker)— = [Kl{Kea) - Kl{Kea)\ 

Ja r I 

+ KeaKQ{Kea)Ki{Kea). (42) 

Inserting these expressions in Equation (I38p 'Dj{k) can be cast as 


PiKl{Kea) 

kf 

1 \ [Jliha) + Jliha)] 



- kiaJQ{kia)Ji{kia)^ 


PeJliha) f 
< \ 

2 [Kli^ea) Kl{Kea)\ 



+ KeaKQ{Kea)Ki{Kea)^. 

(43) 


Furthermore, the amplitudes of improper modes, Af{k), are obtained in a similar fashion by taking 
the scalar product of Equations (fMl) and (l35]) with £,u){r, k). Then we have 











where 


2ujq{ijj) 


( 44 ) 




r±n.\ _ Pi^ h 


2 1,2 ra, 


^.{k) = - 


{ujf ±ig)kirJi{kir) dr 


ki JO 

poo 

PeO^ / {ujf ±ig)ker[CjJi{ker) 

J a 

+ CYYi{ker)] dr. 


(45) 


We can finally present the main expression of this paper. The inverse Fourier transform of 
r, k) with respect to k yields the radial displacement as a function of r and z, 

1 

^r{t,r,z) = —J lrit,r,k)e^^^ dk, (46) 

where ^r{'t,r,k) is given in Equation (I30p . All other perturbations (i.e., magnetic field, velocity, 
and density) are described by analogous formulas. 


5.2. Amplitudes for Specific Initial Conditions 

The expressions presented in Section 15.11 are valid for any axisymmetric initial conditions 
described by the functions f{r,z) and g{r,z) in Equation (f28l) . Here we consider g{r, z) = 0 (i.e., 
the perturbed system is initially at rest) and /(r, z) given by 


fir, z) = ioi’ir) exp(-2;VA^), 


(47) 


with 


V>(r) 


r/a, r < a, 
a/r, r > a. 


(48) 


The parameter A determines the length of the initial tube disturbance. In addition, is the 
maximum radial displacement of the tube boundary at t = 0. In our calculations we consider 
various values of A. 

The Fourier transforms of /(r, z) and g{r, z) with respect to z are 
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fir, k) = ^oV^A'0(r) exp{-A‘^k'^/4:), 
g{r, k) = 0, 

and so the following formulae for Mp{k) and M^{k) can be derived 


(49) 



X 


—a^Qy/TrAujj exp(—A^A;^/4)/sro(«^ea) 

f PiJ2{kia) _ peJoikio) \ 

I k‘f Kj /’ 


(50) 


■^^(k) = -a%V^Awexp(-A^A:^/4)|^^J2(A:ja) 

~k PeiCjJoik^a) + CvYo^k^a)] i. 


(51) 


Therefore, analytical expressions for all the integrals in the amplitudes A^ik) and A^{k) have been 
obtained. Note that if g{r,z) = 0 then A'^{k) = Aj (k) and A^{k) = A~{k). 


6. RESULTS 

The procedure to obtain the perturbations with the help of Equations (l30l) and (H6l) is described 
in Section 5.4 of Paper I and thus is not repeated here. In this work we concentrate on the radial 
displacement at a point on the magnetic tube boundary (r = a). The dimensionless results only 
depend on two parameters, namely the density ratio PijPe (or equivalently the Alfven speed ratio 
UAi/uAe) and A. This parameter determines the length of the initial perturbation, which is of the 
order of 4A; see Equation (j47)l . We use the parameter values pijpe = 4,9 and A = a, 2a. If 
dimensional results are needed, then two more parameters are required: the tube radius, a, and 
Alfv& speed, UAi. 


6.1. Contribution of Proper and Improper Modes 

We hrst consider the radial displacement at the magnetic tube boundary for t = 0. This is 
obtained from Equations (j46p and (I30p . where the contribution of proper and improper modes 
are computed separately. The black line of Figure [3] shows that this numerically computed 
reproduces the Gaussian initial condition of Equations ()28p and ()47p . In addition, Figure [3] also 
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reveals that improper modes (red line) have a large contribution to the initial wave form. For the 
considered parameter values, this contribution is larger than that of proper modes (green line). 
The reason for the large amplitude of improper modes can be found in the wavenumber cut-off of 
the fundamental fast sausage eigenmode, that does not exist for \k\ < kd- Its overtones have even 
larger wavenumber cut-offs. Then, given the initial disturbance of Equation (1471) and its Fourier 
transform of Equation ()49p . that peaks at A: = 0, the contributions to f{r,k) for \k\ < kd cannot 
be accounted for by proper modes. This implies that improper modes must have an important 
contribution to the initial amplitude of the perturbation and, therefore, to the posterior evolution 
of the wave train propagation. Another signihcant feature of Figure [3] is that the contribution 
of both proper and improper modes extends well beyond the length of the tube affected by the 
initial disturbance. The two contributions display spatial variations that added together yield the 
Gaussian shape of Equation (1491) . 

It is worth noticing that improper eigenmodes are more relevant for the present axisymmetric 
disturbance than for the transverse disturbance: compare Eigure [3] and Eigures 4(a) and (b) of 
Paper I. The main difference between fast kink and sausage proper modes is that the fundamental 
proper kink mode exists for all wavenumbers, so that it can account for the main part of the function 
/(r, k). This fact lowers the importance of improper modes in the case of a transverse disturbance 
and also has a second effect: contrary to the case studied here, at t = 0 proper and improper modes 
are confined to the vicinity of the initial perturbation, such as revealed by Figures 4(a) and (b) of 
Paper I. 

From Figure [3] we also see that the contribution of improper modes increases when either A 
is increased or pijPe is decreased. We next explain this behavior. An increase of A corresponds to 
a longer initial perturbation, whose Fourier transform becomes more concentrated around A: = 0. 
Since the fundamental fast sausage mode only exists for |A:| > A:ci, an increase in A results in a 
decrease of /(r, k) for \k\ > A:ci, so that the contribution of proper modes becomes smaller. This is 
corroborated by Figures |D(a) and (b), where the contribution of proper modes to the amplitude of 
at r = a is shown for two values of A. In addition, the role of the density ratio is to modify the 
cut-off wavenumbers, in such a way that a decrease in pijpe leads to a larger wavenumber cut-off 
(see Figure [2|). This then yields a smaller contribution of proper modes: compare Figures IH a) and 
(c). 

Figure 0] also provides information about the relative importance of proper mode overtones. 
Such as expected, the largest amplitude of all proper modes corresponds to the fundamental fast 
sausage mode, whereas its radial overtones possess smaller amplitudes. By comparing panels (a) and 
(b) we see that longer initial perturbations result in a negligible amplitude of all radial overtones. 
We must mention that in Figures [D(a) and (b) only the fundamental mode and its first overtone are 
present (solid and dashed lines), although in panel (b) the first overtone has negligible amplitude, 
so its curve is invisible. Furthermore, in Figure [11(c) the dashed and dotted curves overlap for 
3 < A:a < 4 and may seem to be just one. 
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6.2. Propagation of Wave Trains 

Next, we concentrate in the evolntion of the pulse of Figure[3l^a) . This initial perturbation splits 
in two identical wave trains that propagate symmetrically in both directions along the magnetic 
tube. The early stages of the propagating wave trains are presented in Figure [5] and Animation 1. 
Figures [5]^a)-(c) correspond to the separation of the two wave trains, whereas in Figure [Hd) these 
two features are fully formed. The strong dispersion of the axisymmetric perturbation can be well 
appreciated in this last panel: both wave packets have developed some 6 maxima and minima and 
their length is more than 4 times that of the perturbation at t = 0. Proper and improper modes 
contribute with similar amplitudes to the total signal represented in Figure [5] and Animation 1. A 
secondary effect of wave dispersion is the strong reduction of the wave packets amplitude. In the 
absence of dispersion, the two wave packets of Figure [5|^d) would have amplitude 0.5, whereas their 
actual peak-to-peak amplitude is under 0.3. 

When discussing Figure [3] we explained that, when considered separately, the contributions 
of proper and improper modes affect a section of the magnetic tube much longer than that of the 
initial perturbation. At t = 0 these two contributions simply cancel out in the parts of the magnetic 
tube where the initial excitation is negligible. This situation also takes place for t > 0: although 
proper and improper mode contributions disturb a large part of the magnetic tube, their sum is 
concentrated in a smaller section. Thus, ahead of the propagating wave fronts, where the initial 
perturbation has not yet arrived, these contributions are spatially out of phase and cancel out. This 
is also found at —2 < zja < 2 of Figure [5]^d) because the two wave packets have moved beyond 
this section of the magnetic tube, which is barely disturbed. 

The green and red curves of Figure [5] reveal the separate behavior of proper and improper 
modes. The first ones propagate with group velocity ranging from ~ 0.84uAi (see the minimum of 
the black line in Figure [2Kc)) to UAe = 2uAi. Because of this variation of group velocity the green 
curve displays some dispersion, that becomes more clear for later times (this issue will be discussed 
below with the help of Figure [6|). On the other hand, improper modes propagate at their own 
phase speed and so are subject to phase mixing, that leads to stronger dispersion. This dispersive 
behavior in turn causes a severe amplitude decrease that is clearly visible when comparing the red 
lines of Figures [5]( a) and (d). The asymptotic behavior of improper modes for large t was studied 
in Section 5.3 of Paper I, that deals with a transverse initial perturbation. It was found that 
the contribution of improper modes decays as 1/t. The same applies to the present axisymmetric 
perturbation. 

At t = 20rAi the perturbation of Figure [3K a), whose early development is presented in Figure[5l 
is spread over a long section of the tube of length ~ 30a (compare this value with the initial length 
of the perturbation ~ 4a); see Figured]) a). Here the total perturbation of the radial displacement 
and the contributions of proper and improper modes are represented with the usnal colors (black, 
green, and red, respectively). The vertical lines denote the position of a point initially at z = 0 
and moving at specific velocities: the dash-dotted, dashed, and dotted lines correspond to the 


- 17 - 


maximum group velocity of proper modes (2uAi), the minimum group velocity of the fundamental 
fast sausage mode (0.84uAi), and the minimum group velocity of its first overtone (0.73uAi)- This 
range of group velocities determines the length of the wave train at a given time. The wave train 
shape is controlled by two more features: the dependence Cg{k) of proper modes (black lines in 
Figure[2](c)) and the amplitude of the perturbation (Figure |D(a)). Wavenumbers slightly above the 
cut-off occupy the leading part of the wave packet and have the largest amplitudes, but their Cg 
varies strongly with k. Thus, these wavenumbers are strongly dispersed and contribute moderately 
to the wave train amplitude. On the other hand, wavenumbers ka > 1.6 have group velocities in 
a small range around the minimum of Cg and for this reason the maximum amplitude of the wave 
train is found near the dashed line of Figure[6Ka). A side effect of the wave dispersion is the decrease 
of the wave packet amplitude (compare Figures [5Kd) and[6Ka)). Again, note the cancellation of 
proper and improper mode contributions ahead and behind the propagating wave train. 

We now investigate the effect of the initial disturbance length and the density ratio on the wave 
train dispersion. If the internal to external density ratio is kept fixed and the length of the initial 
disturbance is doubled (Figure [61(b)) the dispersion of the wave train in its constituent frequencies 
does not change because the group velocity of proper modes is the same. Then, the positions of 
maxima and minima in Figures [6Ka) and (b) are the same. Nevertheless, according to Figures [3Ka) 
and (b) and Figures UK a) and (b), a longer initial disturbance corresponds to smaller amplitudes of 
proper modes and for this reason the amplitude of the wave train is much smaller for A = 2a than 
for A = a. 

Finally, we consider the case of a magnetic tube with a larger internal to external density 
contrast {pijpe = 9; Figure [6Kc)). In this case the minimum group velocity of the fundamental 
fast sausage mode and its hrst overtone are similar to those for pijpe = 4 and this implies that 
the position of the vertical dashed and dotted lines in Figures [6Ka) and (c) are similar. But by 
increasing the density ratio from 4 to 9, the maximum group velocity of proper modes increases 
from 2uAi to 3uAi and so the front of the wave train travels considerably faster. Hence, the wave 
packet extends over a much longer section of the magnetic tube, with length ~ 50a. The leading 
part of the wave train has a small amplitude because of the strong gradient of Cg{k) near the cut-off 
wavenumber. In addition, the largest amplitude is found at the back of the wave packet, were, 
again, a wide range of wavenumbers with values of Cg near the minimum group velocity concentrate 
their contribution. 


6.3. Comparison of Signals Caused by Axisymmetric and Transverse Disturbances 

Here we compare the response of the magnetic tube boundary under initial axisymmetric 
and transverse disturbances, respectively described in this paper and in Paper I. To make this 
comparison we consider the wave train passage at a fixed position {z = 20a) and plot the radial 
displacement as a function of time (see Figure [7l in which the left and right columns correspond to 
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axisymmetric and transversqj disturbances). 

Inspection of Figure [7] reveals the very different dispersive properties of fast sausage and kink 
waves. The first ones display rapid oscillations with period ^ 3rAi, whereas the latter have only one 
full oscillation with larger amplitude and period ~ 6rAi. This disparate behavior occurs because 
the fundamental fast kink mode does not have a cut-off wavenumber and its group velocity has a 
more moderate range of variation than that of its sausage counterpart. The length of the initial 
perturbation plays a different role in the propagation of axisymmetric and transverse disturbances. 
In the first case (compare Figures [TKa) and (c)) a longer disturbance imparts less energy to proper 
modes and the wave train disperses more quickly, with a stronger decrease of its amplitude. On 
the other hand, a longer transverse disturbance (compare Figures [T^b) and (d)) imparts less energy 
to the hrst overtone (and so the trailing oscillations are missing in panel (d)) and the fundamental 
mode receives more energy in the long-wavelength regime, for which Cg has a small spread. This 
then results in less dispersion of the wave train. It is also worth noting that the first fast kink mode 
overtone leaves its imprint in Figures E^b) and (f) as small amplitude, short period oscillations 
that arrive at the detection point after the main pulse. Finally, the small amplitude oscillations in 
panels (a) and (c) before and after the wave train passage are numerical artifacts. 


6.4. Wave Leakage 


_ ^The improper modes are cl osely related to leaky quas i-modes studied by, e.g., ICallvl (jl986l . 

2003l h see also the discussion bv lAndries &: GoossensI (j2007l i. The interference of improper modes 


will create a wave front propagating from the tube in the radial direc tion. The leaky quasi-modes 
describe the so-called intermediate asymptotic of this wave front (see iBuderman fc RobertsI 120061 . 
for a discussion of the role of leaky quasi-modes in the case of kink waves). 


Wave radiation after the initial excitation can be well appreciated by plotting snapshots of the 
dependence of with respect to r for several times (Figure [8|). Panel (a) shows the imposed radial 
prohle of the radial displacement at t = 0, given by Equation (H8|) . with improper modes (red curve) 
having a more important contribution than proper modes (green curve). The subsequent temporal 
evolution of shows oscillations inside the magnetic cylinder (0 < r/a < I) and in its vicinity, 
together with a wave pulse propagating outwards in the radial direction, that is specially visible in 
panel (d) as a depression in in the range 3 < r/a < 7. This propagating feature is clearly caused 
by improper modes because the amplitude of proper modes at its position is negligible. 

The wave emission process has been presented here by doing the cuts of Figure [8] at the center 


^This case was studied in Paper I, although the external variation of '(/’(r), which gives the radial profile of the 
initial perturbation, was different from that used in the present work (Equation (08}). In order to make a comparison 
with the results of Paper I, the proper and improper mode amplitudes of transverse disturbances are presented in 
Appendix o for the same ip{r). 
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of the initial perturbation (z = 0), but as far as the amplitude of improper modes is non zero, it also 
takes place at other positions along the magnetic tube when they are affected by the propagating 
disturbance. 


7. DISCUSSION 


In this paper we have used a method based on Fourier integrals to solve the initial value problem 
of impulsively excited axisymmetric (i.e., sausage) perturbations in a uniform, dense, cylindrical 
magnetic tube. The linear and /3 = 0 approximations have been assumed. This method gives 
analytical expressions for the perturbations (plasma displacement, density, velocity, and magnetic 
field) that require the numerical computation of the two integrals in Equations (1301) and (14611 . The 
accuracy of the results depends exclusively on the precision of the numerical computation of these 
integrals. An advantage of this method is that one can select which perturbations will be computed 
and at which spatial position(s) and time(s). 

The impulsive generation of fast sausage wave trains in a dense magnetic tube has also been 
studied before by numerically integrating the line arized MHD equations. Mos t of these works use 


a Ca r tesian slab to represent the magnetic tube (IMurawski fc Roberts 


20041 . 12 OO 5 I : Ijelmek &: Karlickv 2010|; iNakariakov et al 


2012 


1993a 


^1 


Meszarosova et al. 


Nakariakov et al. 


2 OI 4 II . In the solar 


context, only ISelwa et al.l (j2004l ) have considered this problem in cylindrical geometry, although 
they did not take into account the dispersive behavior of the propagating wave packet. 

One of the key features of the present paper is the importance of the group velocity in the 
dispersion of a c oncentrated perturba tion. Concerning fast sausage wave trains, this was already 
acknowledged by 


Nakariakov et al. 


Roberts et al. (1984) and later by other authors (e.g., Nakariakov fc Robertslll995 


20041 1 . In addition to this, we have emphasized the role of the length of the initial 


perturbation and the internal to external density ratio in the wave packet dispersion (Figures [6] 
and ED. We have found that longer initial perturbations have much lower amplitudes because the 
initial excitation imparts less energy to proper modes. In addition, larger values of the density 
ratio lead to a wider wave packet that nevertheless has a larger amplitude because most of the 
oscillatory energy is concentrat ed in its trailing part. S ome of these results have also been found 
in a non-uniform dense slab bv INakariakov et al.l (120051 ) 


We are now concerned with the observational implications of our work. One of the main con¬ 
clusions of Paper I is that a spatially concentrated, transverse disturbance of a magnetic tube gives 
rise to non-negligible density perturbations that take their maximum value at the tube bound¬ 
ary. The reason for these density perturbations is that fast kink waves are compressive for lon¬ 
gitudinal wavelengths of the order of the tube radius. The axisymmetric wave trains studied in 
the present work also generate density perturbations caused by the compressive character of fast 
sausage waves. To de t ermin e whether the wave trains propagating along coronal loops studied by 
Williams et al.l (120011 . [2002); iKatsivannis et al.l (120031 ) (and similar events detected in solar radio 
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bursts, see iMfezarosova et alJl2009l . I 2 OIII : iKarlickv et alJl2013l ) are caused by fast sausage or kink 
excitations it is necessary to carry out a forward modelling of the magnetohydrodyna mic perturba¬ 


tions to determine how they would appear when observed with available instruments. ICooper et al 


(|2003al jbf) assumed that the observed emission intensity in optically thin coronal lines is propor¬ 
tional to the line integral along the line-of-sight (LOS) of the electron number density squared. 
They considered fast sausage and kink waves with fixed wavelength and different orientations of 
the magnetic tube axis with respect to the LOS. They concluded that the two wave modes generate 
similar amplitudes of integrated density variations. Nevertheles s, the integrated intensity alone 
does not govern the line intensity and formation of emission lines. iDe Moortel &: Bradshaw! (j2008l l 
took into account that it is actually governed by a complex, non-linear interaction between the 
density squared, the ionisation balance and emissivity, and the instrument response function. The 
important point here is that the ionisation equilibrium timescale of Fe X and Fe XII ions in a coronal 
environment is of the order of a few to tens of seconds. The effect of the ionisation balance can be 
ignored only if the wave period is much longer than this timescale. Otherwise, observed intensity 
perturbations do not necessarily behave like the theoretical temperature a nd density. Therefore , 
to m a ke a theoretical in t erpre tation of the ~ 6 s oscillations described by IWilliams et al.l ([200l|, 


2002 ! ) ; iKatsivannis et al.l (1200.31) a careful forward modelling of imp ulsively e xcited fast sausage and 


kink perturbations is required. lAntolin fc Van Doorsselaerel (l201,3l ) (see also IR eznikova et al.ll2014l ) 
studied the influence of LOS angle and spatial, temporal, and spectral resolution of observing in¬ 
struments on fast sausage wave observables. Under the assumption of equilibrium ionisation, these 
authors concluded that the intensity variation can be very low (< 4% for Fe IX 171 A) or significant 
(< 35% for Fe XII 193 A). Most of the intensity variation, however, disappears when considering 
the radiative relaxation times of the ions. Another important conclusion of this work is that the 
intensity modulation has strong dependence on the LOS angle and on the spatial and temp oral 
resolutions. Our future aim is to repeat the calculations of lAntolin fc Van Doorsselaerel (j2013l ) for 
impulsive initial excitations of the fast sausage and kink types. 
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A. FREQUENCY OF PROPER MODES FOR LARGE k 

Our purpose here is to study the solutions to the dispersion relation 

Ji{kia) ^ Ki{Kea) 


kiaJQ{kia) KeaKQ{Kf.a) 


(Al) 


in the limit of very large longit udinal wavenumber, i.e., for fca 3> 1. We make use of the asymptotic 
formulae for large arguments (jAbramowitz fc Stegunlll964l ) 


(A2) 


Ko{x) ~ Ki{x) ~ 


Figure [2] shows that, for ka 3> 1, oj/k is not close to UAe and so KgO S> 1. Then, the second term 
in the dispersion relation is 

« 1. (A3) 

KeaKo[Kea) KgO 

Again using Figure [2] we see that, for ka 3> 1, co/k is close to UAi and so ki is not large. Since 
the second term in Eqnation (lAip is small, the first one also has to be small. This is only possible 
when kia is close to one of the zeros of the function Ji(x), which are jij, j = 1, 2 ,... As a resnlt, 
we have 


= fc^UAi I 1 + 




(A4) 


Y Pek'^a? j 

From this expression, the phase speed and group velocity for large longitudinal wavenumber are 


Cph UAi 1 + 


\ 

2pek‘^a?‘ j ’ 


1 - 


_Pi^ 

2pek‘^a? 


(A5) 


We conclude that as ka tends to infinity, Cph tends to vai from above and Cg tends to vai from 
below. 


B. CALCULATION OF q{u:) 


Using Equation (j26[) we write the expression for the scalar product of two improper eigenfunc¬ 
tions as 

Pia^k^k'J^ 


where 


/ Ji{kir)Ji{k'ir)rdr 

Jo 

poo 

+ Pea'^kek'^ / Fi{ker)Fi{k'^r) r dr, 
J a 

F,{x) = CjMx) + CyY,{x). 


(Bl) 

(B2) 
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Now we use the formulae 



where ji{z) and yi{z) decay as 2 ^ when z ^ 00 . 
rewrite Equation ()B1D as 


Itt) [l+ji{z)] , 

In) [1 + yi{z)] , (B3) 

Using these formulae, after some algebra, we 


{L,L') = T{uj,uj') + T(a;,a;'), 


(B4) 


where 


r(w, — El (w, + r2 (w, + r3(u;, w^), 

„.„4L2jL/2 pa 

ri{u],u]') = Ji{kir)Ji{k[r)rdr, 

Jo 


(B5) 

(B6) 


T2{uj,u;') = 


Pea^^/KK 


TT 


poo 

/ {[C'jCSii(A;'r) + CjC;.yi(A:;r)] 

J a 


X cos {kef — 1^) 

+ [CYC'jh{k'er) + CYC'yyi{k'eT)\ 
X sin {kef — Itt) 

+ [CjC'jhiker) + CYC'jyiiker)] 
X cos {k'^r — Itt) 

+ [CjC{. hiker) + CYC{.yhker)] 
X sin {k'eV — Itt) } dr, 


(B7) 


T3{uj,uj') 


PeO^ 

n^/kek'e 


poo 

/ [CjCjhiker)jiiker) 

J a 


+ CjCYh{ker)yi{k'er) 

+ CYC'jyi{ker)ji{ker) 

+ CYCYyi{ker)yi{ker)] dr, 


T{UJ,UJ') = ^^^y/kek'e 
7T 
poo 

X / {{CjC'j + CYC{.)cOs[r{ke-k’e)] 

~ {CjC'y + CyC'j) cos \r{ke + fcg)] 

- {CjC'y - CyC'j) sin [r{ke - k'h] 

+ {CjC'j - CyC'y) sin [r{ke + k'e)] } dr. 


(B8) 


(B9) 
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As it is explained in Paper I, r(a;,a;') is a regular function of u and a;'. Repeating the same 
calculation as in Paper I we obtain that 

(4,0) = Q(a;,a;') + - u'), (BIO) 

OJ 

where Q{uj,uj') is a regular function and 

W{00,00') = {CjC'j + CyC'y) cos [a{K - K)] ■ (Bll) 

Using the same arguments as in Paper I we prove that Q{oo,oo') = 0 and W{oo,oo')5{oo — 00') = 
W{oo,oo)6{oo — 00'). Then we obtain 

,(„) = (c3 + cj) ^ (B12) 

OJ 

C. AMPLITUDES OF FAST KINK MODES FOR SPECIFIC INITIAL 

CONDITIONS 

In Paper I the evolution of an initial transverse displacement of the magnetic tube was studied. 
The radial variation of the initial disturbance outside the magnetic cylinder was taken as exp(—(r— 
with I a prescribed length. In the present paper this external radial variation is a/r. We 
now derive expressions for the amplitudes A^{k) and A^{k) for a transverse perturbation whose 
initial radial distribution is for r > a. Note that for this type of perturbation we still 

use = 1 for r < a. 

We start with the proper mode amplitudes A^(k), that are given by Equation (fMl) with 


Afj{k) = a^o^/^Aooj exp{—A‘^k‘^/4:)Ji{kia) 

Pj_ 

k? 


^ j lY ( \ P^ 

X < —Ai(Kea) - 


TT 


- - + K,eaKo{Kea) 


+ ^K{Kea) 


(Cl) 


where 


vr 


^k{x) = -X [Ko{x)Li{x) + Ki{x)Lo{x)] . 


(C2) 


This function is introduced by Rosenheinrichl (j2014l '). p. 6, and appears when obtaining the primitive 
of Ki{x)/x, see IR osenheinrichl (120141 ). p. 19. Furthermore, Lo(x) and Li(a:) are the modified Struve 
functions. 
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The denominator of Equation ()36l) contains 


+ {kfa^ - 2 ) jf{kia)} 

+ 


2At4 


+ ~k 2) (/te®)} ■ 


(C3) 


To derive this result the primitive of xK'i{x) is required. We have used Equation (l39]l to express 
K[{x) in ter ms of Kn(x) and K^( x). Next, the primitive has been calculated with the help of 
expressions in iRosenheinrichl (1201411 . pp. 143 and 160, with Io{x) and Ii{x) replaced by Kq{x) and 
Ki{x). 

Finally, Equation (1441) gives the amplitudes A^{k), where the numerator is 


Afuj{k) = exp(—A^/c^/4) 


/ Pikjo 

I kf 


+ Pea 


-Ji{kia) 

Cj{l - keaJo{kea) - ^{keo)) 

+ Cy{- keaYo{kea) - ^yikea)) 


(C4) 


where the primitive of Ji{x)/x is given in IRosenheinrichl ((2014), p. 19, in terms of the function 


TT 


4>(x) = -X [Ji(x)Ho(x) - Jo(x)Hi(x)] , 


(C5) 


with Ho(x) and Hi(x) the Struve functions. The primitive of Y\{x)/x is derived in the same way 
by substituting Jo{x) and Ji(3:) by To( 2 ;) and Yi{x) and 4>(a;) by 4>y(x). This function is defined 


as 


TT 


$y(x) = -X [yi(x)Ho(x) - yo(x)Hi(x)] . 


(C6) 


The numerical evaluation of the Struve and modified Struve functions has been done with the 


software developed bv iMacLeodI (jl996l ). 
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(a) 





Fig. 2.— (a) Frequency, (b) phase speed, and (c) group velocity vs. the longitudinal wavenumber for 
the fundamental fast sausage mode (solid line) and its first overtones in the radial direction (dashed, 
dotted, and dash-dotted lines). Black and red colors respectively correspond to the equilibrium 
density ratio PijPe = 4 and 9, i.e., to the Alfven speed ratio UAe/'i’Ai = 2 and 3. The purple lines 
in (a) and (b) correspond to the frequency u = kv^e- To obtain the dimensionless quantities, the 
magnetic cylinder radius (a), Alfven speed (uAi), and transit time (tai = a/uAi) have been used. 
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z /a 

Fig. 3.— Initial shape of the magnetic tube boundary, given by Equations (H7|) and (H8|) with r = a. 
The radial displacement of the tube boundary is plotted at f = 0 as a function of the distance along 
the magnetic tube (black line). The contributions of proper and improper modes are shown as 
green and red lines, respectively, (a) Density ratio PijPe = 4 and A = a; (b) pi/pe = 4 and A = 2a; 
(b) Pij Pe = 9 and A = a. According to Equation (|T7ll the length of the initial disturbance is of the 
order of 4A. 
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Fig. 4.— Proper mode contribution to the radial displacement as a function of the longitudinal 
wavenumber. When inserted in Equation (j30p . the product A^{k)^j{a,k) determines the contribu¬ 
tion of each proper mode to the radial displacement for a fixed k. (a) PijPe = ^ and A = o; (b) 
Pi!Pe = 4 and A = 2a; (c) pij Pe = 9 and A = o. Solid, dashed, and dotted lines correspond to 
the fundamental mode and its first two overtones. The vertical dotted lines mark the wavenumber 
cut-offs. 
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Fig. 5.— Shape of the magnetic tube boundary for several times (shown at the top of each frame). 
The radial displacement at r = a is plotted as a function of z. This figure shows the early evolution 
of the initial pulse of Figure[3](a). Black lines correspond to while green and red lines correspond 
to the contributions of proper and improper modes, respectively. The density ratio is pijPe = 4 
and the length of the initial disturbance is of the order of 4a. 

(An animation of this figure is available.) 
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Fig. 6.— Shape of the magnetic tube boundary at t = 20rAi. The three panels correspond to the 
wave trains caused by the initial disturbances of Figure [3l with the same meaning of line colors. The 
vertical lines give the position of a point initially at z = 0 and traveling at a constant speed equal 
to the maximum group velocity of proper modes (dash-dotted line), the minimum group velocity 
of the fundamental proper mode (dashed line), and the minimum group velocity of the first proper 
mode overtone (dotted line). 
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(a) (b) 








Fig. 7.— Radial displacement of the magnetic tube boundary at a distance z = 20a from the initial 
perturbation. Left column: axisymmetric (sausage) perturbation; right column: transverse (kink) 
perturbation, (a) and (b) = 4, A = a; (c) and (d) = 4, A = 2a; (e) and (f) pi/pe = 9, 

A = a. The vertical lines have the same meaning as in Figure [H 
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Fig. 8.— Radial displacement of the plasma as a function of radial distance at z = 0 and for 
different times, shown above each frame. Black lines correspond to the displacement while 
green and red lines correspond to the contributions of proper and improper modes, respectively. 
The density ratio is pi/pe = 4: and the length of the initial disturbance is of the order of 4a. 




















